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При анализе надежности конструкций из полимерных и полимерных композитных материалов применяется вероятностный под- 
ход, поскольку для реальных материалов всегда существует некоторый разброс количественных характеристиках их свойств, 
для конструкций — отклонения размеров от их номинальных значений, для нагрузок — отклонения от средних эксплуатационных 
значений. Обработка массива параметров напряженно-деформированного состояния конструкций, полученных в результате 
численных и натурных экспериментов, проводится с помощью методов теории вероятности и математической статистики. 


Введение 

При анализе параметров напряженно-деформи- 
рованного состояния (НДС) конструкций и оценке 
их работоспособности распространенным является 
подход, который можно назвать детерминистским. 
Принимается, что все параметры, которыми опре- 
деляется НДС конструкции, являются определен- 
ными с известной точностью величинами. Для за- 
дач анализа НДС деталей и конструкций это три 
группы параметров. Первая из них определяет 
свойства материала, вторая - геометрию конструк- 
ции, третья - способ приложения и интенсивность 
приложенных нагрузок. Все эти характеристики в 
той или иной степени носят случайный характер, а 
конкретные их величины, задаваемые в детермини- 
стском подходе, являются некоторыми средними 
приближенными значениями. Существующий 
разброс параметров и их отклонение от средних 
значений учитывается введением коэффициента 
запаса прочности, чем компенсируется неопреде - 
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ленность информации о точных значениях. Можно 
говорить, что величина коэффициента запаса проч- 
ности - это характеристика уровня знания (вернее, 
незнания) точных значений этих параметров. 

Возможность разрушения изделия, таким обра- 
зом, носит вероятностный характер, и она должна 
оцениваться соответствующими количественными 
характеристиками. Вероятностный характер имеют 
не только выше перечисленные параметры матери- 
ала, конструкции и нагрузки, но и сами критерии 
разрушения, поскольку они представляют собой 
обработку некоторого массива экспериментальных 
данных. При оценке вероятности безотказной ра- 
боты (ВБР) необходимо сопоставлять расчетные 
или экспериментальные данные о НДС конструк- 
ции во всех ее точках, имеющие вероятностный ха- 
рактер, с критериями прочности, имеющими такой 
же смысл. Подобные подходы разрабатываются, 
например, в [1, 2]. Отличительные особенности 
предлагаемой работы заключаются в учете разбро- 
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са параметров материала не только от опыта к опы- 
ту, но и по самим конструкциям, что связано с тех- 
нологическими особенностями их изготовления. 

Физическая и математическая постановка 

Качество изделий из полимерных композитных 
материалов (ПКМ), получаемых различными техно- 
логическими способами, определяется свойствами 
его фаз, зависящими как от рецептуры, так и от ре- 
жимов переработки материала. В большинстве слу- 
чаев операции изготовления ПКМ и изделия из не- 
го совмещаются при использовании различных тех- 
нологий изготовления изделий: свободным литьем 
и литьем под давлением, прессованием, экструзией 
и т.д. Нестабильность параметров технологических 
процессов приводит к тому, что деформационно- 
прочностные свойства изделий и их геометрические 
размеры являются не детерминированными, а слу- 
чайными значениями. Более того, в разных точках 
по объему пресс-формы отличаются уровни давле- 
ния, температурный режим и т.д. Случайными вели- 
чинами в условиях эксплуатации являются и внеш- 
ние воздействия, например, снеговая или ветровая 
нагрузка, скачки давления в трубопроводе при сра- 
батывании запорной арматуры или включении и 
выключении насосов и т.д. Для описания поведения 
конструкций всем прочностным и геометрическим 
характеристикам, а также внешним воздействиям 
придается вероятностный характер. 

Для построения полей ВБР конструкций из 
ПКМ с учетом нестабильности их свойств опреде- 
ляются параметры НДС конструкции, которые но- 
сят вероятностный характер, строятся кривые 
распределения значений эквивалентных напряже- 
ний и/или деформаций в каждой точке конструк- 
ции и сравниваются с экспериментально получен- 
ными предельными напряжениями, имеющими 
также вероятностный характер. 



Рис. 1. Варианты распределений эквивалентных деформаций 
и предельных значений в точке конструкции е ш 

Рассмотрим возникающие в связи с этим вариан- 
ты (рис. 1). Считаем для конкретности, что критери- 
ем, по которому можно судить о работоспособности 
материала, является так называемая эквивалентная 
деформация. На оси эквивалентных деформаций е ш 
отрезком АВ обозначим область, в которую попада- 
ют все эквивалентные деформации в данной точке 
конструкции при натурном или численном модели- 


ровании НДС изделия. При численном анализе это 
означает проведение ряда расчетов, в которых часть 
параметров или все параметры, определяющие НДС 
конструкции, являются случайными величинами. В 
натурном эксперименте это означает проведение ря- 
да испытаний одного и того же или нескольких по- 
добных изделий. На этой же оси отрезком АД обоз- 
начим интервал эквивалентных деформаций е ж *, со- 
ответствующих разрушению, или пределу текучести 
материала при экспериментальных исследованиях. 

Тогда в случае а (рис. 1) в данной точке 
конструкции разрушение исключено, в случае Ь 
разрушение произойдет обязательно. В случае с об- 
ласть "перекрытия", отрезок А, В, есть не что иное, 
как характеристика (вероятность) разрушения 
(точнее, она определяется отношением А,В/АВ), а 
в случае сі - область АВ, (АВ,/АВ) - вероятность не- 
разрушения. 

Эти оценки справедливы, если вероятность по- 
падания эквивалентных деформаций в данной точ- 
ке конструкции в любое место внутри отрезка АВ 
одинакова, а уровни е ж *, отвечающие разрушению 
материала при его испытаниях, равномерно расп- 
ределены по отрезку А, В!. Экспериментально уста- 
новлено, что на самом деле распределения этих па- 
раметров подчиняются нормальному закону. В 
этом случае для вариантов а и Ь (рис. 1) все рассуж- 
дения сохраняют свою силу, а для вариантов с и сі 
нужно оценивать так называемую область "перек- 
рытия" кривых распределения е жв и е„*. 

Область "перекрытия" (рис. 2), полученная в ре- 
зультате сопоставления кривых распределения эк- 
вивалентных деформаций и их предельных значе- 
ний, и является областью вероятного отказа рабо- 
ты конструкции [2] . 

Де экв ),Ле* экв) 



Рис. 2. Пример распределений: 1) предельной и 2) реальной 
эквивалентных деформаций 

В работе предлагаемый подход иллюстрируется 
на примере построения поля вероятностей безот- 
казной работы для конструкции в виде участка тру- 
бы, выполненной из полимерного материала на ос- 
нове полипропилена, используемого для производ- 
ства трубопроводов различного технологического 
значения - транспортировки воды, газа и т.д. Од- 
ним из значимых параметров, определяющим НДС 
изделия, является модуль упругости материала. 
Экспериментально установлено, что от наблюде- 
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ния к наблюдению значения модуля упругости от- 
личаются от его среднего значения в пределах 20 %. 
При математическом моделировании случайное 
изменение модуля упругости от наблюдения к наб- 
людению задается с помощью датчика псевдослу- 
чайных чисел [3]. Согласно экспериментальным 
данным, это изменение подчиняется нормальному 
закону распределения. В качестве эквивалентной 
деформации, по которой можно судить о работос- 
пособности материала в конструкции, принимает- 
ся интенсивность деформаций. В силу осевой сим- 
метрии задачи рассматривается расчетная область, 
представляющая собой часть осевого сечения тру- 
бы, и решается осесимметричная задача теории уп- 
ругости (рис. 3). 
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Рис. 3. Расчетная схема нагруженного участка трубы 


В качестве внешних нагрузок принимаются од- 
новременно приложенное внутреннее давление Р, 
и осевое сжатие Р,. Для упругой трубы под действи- 
ем внутреннего давления известно решение Ламе 
[4], которое в данном случае служит для тестирова- 
ния программы. 

В каждой точке контура области АВСІ.) необходи- 
мо ставить по два условия (задача является двумерной 
по пространственным переменным), задаваясь век- 
тором напряжений, вектором перемещений или по 
одной из разнонаправленных компонент этих векто- 
ров. В примере использованы следующие условия. 

На линии АВ ставятся условия скольжения 
вдоль жесткой стенки: 

и = 0, сг г = О, 

где и - радиальное смещение, <т г - касательные 
напряжения. 

На линии ВС напряжения отсутствуют: 
сг„ = <т г = О, 

где ст„ - нормальные напряжения. 

На линии СБ: ст„ = -Р 2 , сг г = 0. 

На линии АБ: ст„ = -Р,, <т г = 0. 


Метод решения 

Численная реализация задачи определения па- 
раметров НДС конструкции для каждого конкрет- 
ного варианта распределения его свойств по объему 


изделия проводится методом конечных элементов с 
использованием разбиения расчетной области на 
конечные элементы треугольной формы. Метод ос- 
нован на использовании соотношения: 

ягт м 

= [кП{Щ + {Г е) }, (1) 

д{Щ 

где ГГ 11 - потенциальная энергия отдельного эле- 
мента системы, {Ц} - узловые перемещения, 
{7 м } - вектор нагрузки, [СД - матрица жесткости 
отдельного элемента, представляющая собой объ- 
емный интеграл вида: 

[к (е) ]= [ [5 (е) ] 7 '[22 (е) ][5 (е) ]г/С. (2) 

ѵ (е) 


Для разбиения расчетной области используют- 
ся треугольные конечные элементы с шестью ком- 
понентами узловых перемещений. Для уменьше- 
ния вычислительных ошибок треугольные элемен- 
ты выбираются возможно более правильной гео- 
метрической формы, без тупых углов. 

Уравнения для элемента записываются в цили- 
ндрической системе координат. Компоненты пере- 
мещения и, ѵ аппроксимируются внутри треуголь- 
ного элемента соотношениями: 
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Соотношение (1) содержит три функции фор- 
мы, которые в линейном случае имеют вид: 


М і=—М + Ь і г + с і 2), 
2 А 


где А - площадь треугольного элемента, с,- - 
коэффициенты, вычисленные по значениям узло- 
вых координат. 

Для осесимметричной задачи используются ге- 
ометрические соотношения: 


ди дѵ и ди дѵ 

— 5 ~ ~ я 7Г~ • 

дг дг г дг дг 


(5) 


Матрица упругих характеристик [Д] в случае 
осесимметричной задачи и для изотропного мате- 
риала имеет порядок 4x4: 
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Дифференцируя (3) и используя выражения для 
функций формы, получим соотношения для де- 
формаций: 

{8} = [В]-{Щ, (6) 

где [5] - матрица, которая содержит коэффициен- 
ты, являющиеся функциями координат. Матрица 
жесткости [В\ для элемента представляет собой 
объемный интеграл и содержит коэффициенты, 
которые являются функциями координат г и г, и 
поэтому не могут быть вынесены за знак интеграла. 
В [5] предлагается использовать следующую фор- 
мулу для матрицы жесткости: 

[к м ] = [В М ] т ■ [Г> (е > ] ■ \В іе ’ ] • 2 п ■ г А , (7) 

где [Б 1 ' 1 ] - матрица коэффициентов, вычисленная 
по значениям г и % в центре элемента для строки, 
содержащей функции формы. Матричное уравне- 
ние для ансамбля элементов имеет вид: 

= ( 8 ) 

где [ К\ - матрица жесткости, которая собирается из 
матриц жесткости элементов, {/} - глобальный 
вектор-столбец нагрузки. 

е=\ е=\ 

Соотношение (8) представляет собой матрич- 
ную форму записи системы линейных алгебраичес- 
ких уравнений относительно компонент вектора 
перемещений во всех узлах конечно-элементной 
сетки. Эта система линейных уравнений решается 
методом Гаусса. Затем деформации во всех ячейках 
сетки определяются по соотношению (6), а напря- 
жения в элементах вычисляются по закону Гука 
{<т}=[Д]{е}. 

Определение вероятностей безотказной работы 

конструкции 

Задача вероятностного расчета сводится к мно- 
гократному определению параметров НДС, кото- 
рые являются случайными величинами. В результа- 
те обработки полученного массива данных опреде- 
ляются характеристики параметров НДС - матема- 
тическое ожидание и среднеквадратичное откло- 
нение оу Они оцениваются для значений интенсив- 
ности деформации в каждой конечно-элементной 
ячейке расчетной сетки по соответствующим фор- 
мулам теории вероятностей, что позволяет постро- 
ить функции плотности распределения [2]. 

Плотность нормального распределения интен- 
сивности деформаций 5 имеет вид 

і Г і С VI 

ч 1 I 5- /Л, 

/» = /= ех Р -- ,-®<5<оо. 

а, У 2л 2^ <т 5 ) 

Аналогичному закону подчиняется плотность нор- 
мального распределения предельных деформаций А 


1 Г \( 8 -и V" 

/, (В) = ^^ехр — — , -со < 5 1 < да, 

где /и, - математическое ожидание интенсивности 
деформаций, <т 5 - среднее квадратичное отклонение 
интенсивности деформаций, д, - математическое 
ожидание предельных деформаций, - среднее 
квадратичное отклонение предельных деформаций. 

Введем случайную величину у=5-я. Принима- 
ем, что случайная величина у имеет нормальное 
распределение с математическим ожиданием 

и средним квадратичным отклонением 



Тогда ВБР можно выразить через у как 


К = Р(у > 0 ) 




сіу. 


Если х=(у-ц у )/а у , то йх=йу. При у = 0 нижний 
предел случайной величины х имеет вид 

х = 

а при у— >+со верхний предел х->со. Следовательно, 






(9) 


Ясно, что х=(у-/л у )/<т у является нормированной 
случайной величиной, распределенной по нор- 
мальному закону. Следовательно, вероятность бе- 
зотказной работы можно найти с помощью таблиц 
функции нормального распределения. 

Формулу (9) можно представить в виде: 



где Ф - интеграл Лапласа. 

Полученная функция плотности распределения 
вероятностей значений интенсивности деформа- 
ций е і подчиняется нормальному распределению. 
Функция плотности распределения вероятностей 
для экспериментально полученных предельных 
значений интенсивности деформаций е* имеет 
вид, показанный на рис. 2 (кривая 1). Область пе- 
ресечения этих кривых является областью отказа. 
Величина вероятности отказа (безотказной рабо- 
ты) вычисляется по приведенной формуле (10). 
Определив значения х, можно с помощью таблицы 
нормального нормированного распределения оп- 
ределить величину вероятности безотказной рабо- 
ты конструкции. 
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Результаты расчетов 

Сопоставляя полученные деформации с пре- 
дельными экспериментальными значениями, име- 
ющими также случайное распределение (с характе- 
ристиками Д5=0,02, <75=0,001 для ПКМ на основе 
полипропилена), получаем в каждой точке трубы 
по ее толщине значение функции К (ВБР) матери- 
ала. Из анализа решения следует, что наибольшие 
деформации получаются на внутренней стенке тру- 
бы, а вероятность безотказной работы ее в этих 
точках является наименьшей (рис. 4). Если в какой 
либо точке конструкций вероятность безотказной 
работы меньше требуемого (нормативного) значе- 
ния, то можно говорить о вероятности разрушения 
конструкции в целом. В качестве нормативного 
значения вероятности безотказной работы можно 
принять, например, характеристику, связанную с 
отношением стоимости планового ремонта изде- 
лия к стоимости устранения последствий аварий- 
ного разрушения [6]. 

Анализ влияния трех параметров (толщины 
стенки трубы г, внутреннего давления Р х и осевого 
сжатия Р 2 ) на напряженно-деформированное сос- 
тояние конструкции и соответственно на вероят- 
ность безотказной работы К показывает, что наибо- 
лее значимым параметром является А г, далее по 
убывающей Р х и Р 2 . Это означает, что при произво- 
дстве полимерных труб их толщина должна конт- 
ролироваться наиболее тщательно. 
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Рис. 4. Распределение математического ожидания интен- 
сивности деформаций /а, (1, 2) и вероятности безот- 
казной работы Р (3, 4) по толщине стенки трубы г; 
кривые 1, 4) соответствуют уровню внутреннего дав- 
ления Р, = 4,5 МПа и осевого сжатия Р 2 = 1 МПа, кри- 
вые 2, 3)Р=5 МПа, Р г =2 МПа 

Таким образом, ответ на вопрос о надежности, 
или о возможном разрушении конструкции, при 
рассматриваемом подходе получается в виде вероят- 
ности разрушения (неразрушения) материала во всех 
точках изделия, что можно представить в виде полей 
распределения ВБР по всему объему конструкции. 
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